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FOREWORD 


For  the  past  twenty-four  years,  the  Naval  Surface  Weapons  Center  has  been 
involved  in  a  research  and  consulting  effort  for  both  the  Energy  Research  and 
Development  Administration  and  the  Nuclear  Regulatory  Commission  concerned 
with  the  study  of  reactor  vessel  response  to  hypothetical  core  accidents  and 
other  types  of  dynamic  loading  events.  As  part  of  this  effort,  analyses  of  jet 
forces  produced  by  ruptured  steam  process  pipes  on  neighboring  walls  of  reactor 
buildings  have  been  performed.  This  report  presents  the  computed  results  for 
impact  loads  from  the  ruptured  process  pipe  on  the  inner  wall  of  a  surrounding 
concentric  guard  pipe  for  geometries  consistent  with  reactor  plant  design 
guide! ines. 
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CHAPTrR  1 
INTRODUCTION 


1.1  BACKGROUND 

Since  1956  the  Naval  Surface  Weapons  Center  has  been  involved  in  a  research 
and  consulting  effort  for  both  the  Energy  Research  and  Development  Administration 
and  the  Nuclear  Regulatory  Commission  (NRC)  concerned  with  the  study  of  reactor 
vessel  response  to  hypothetical  core  accidents  and  other  types  of  dynamic 
loading  events.  In  an  earlier  effort,  NSWC  provided  to  NRC  preliminary  calcula¬ 
tions  of  pressure  loads  produced  on  neighboring  walls  by  steam  pipe  ruptures. 

The  task  reported  herein  is  a  continuation  of  these  calculations  for  the 
purpose  of  determining  the  pressure  loads  very  close-in  to  the  site  of  steam 
pipe  rupture.  Specifically,  a  series  of  hydrocode  calculations  were  done  to 
determine  the  fluid  flow/pressure  fields  in  regions  bounded  by  an  inner  steam 
pipe  and  an  outer  guard  pipe  following  simulated  ruptures  of  the  high-pressure 
inner  steam  pipe. 

1.2  PROBLEM  STATEMENT 

A  process  pipe,  filled  with  saturated  steam  at  8.3  MPa  (1200  psia),  was 
assumed  to  rupture  instantaneously  in  either  the  circumferential  or  longitudinal 
direction  (See  Figure  1).  The  steam  flowed  out  from  the  rupture,  driving 
ambient  air  ahead  of  it,  and  impacted  the  inside  surface  of  a  guard  pipe  which 
surrounded  the  process  pipe.  The  desired  result  was  the  pressure  distribution 
along  the  guard  pipe. 

The  analysis  was  performed  in  two  phases.  In  the  first  phase,  the  proper 
equation  of  state  for  high-pressure  saturated  steam  was  evaluated,  the  appropriate 
finite  difference  grid  scale  for  the  problem  geometry  and  flow  conditions  was 
determined,  and  a  particular  hydrocode  for  performing  the  calculations  was 
selected.  The  second  phase  of  the  task  determined  the  time-dependent  internal 
pressure  distribution  on  the,  guard  pipe  as  a  function  of  (a)  direction  of  the 
pipe  rupture  (circumferential  or  longitudinal),  (b)  width  of  the  pipe  rupture, 
and  (c)  separation  distance  or  gap  between  the  process  and  guard  pipes.  The 
parameters  for  the  seventeen  process/guard  pipe  configurations  that  were 
computed  are  listed  in  Table  i. 
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FIGURE  1.  GEOMETRY  OF  BREAKS 
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CHAPTER  2 

ANALYTICAL  SOLUTIONS 


For  a  free  compressible  jet  expanding  to  ambient  pressure,  the  total  jet 
thrust  for  an  ideal  gas  as  given  by  Reference  1  (slightly  modified)  is 


T/A„  -  C.PA  -  P 
e  10  * 

where  T  »  Jet  thrust  =  Total  force  on  a  large  normally-impacted  plate 
Ae  »  Effective  area  of  process  pipe  break 
Pq  *  Pipe  reservoir  pressure 
P^  *=  Ambient  pressure 
C,  -  (1  +  y)(2/(y+  1))y/(y'1) 
y  -  Specific  heat  ratio 

The  coefficient  C,  is  a  weak  function  of  v»  being  1.255  for  y  ~  1.3  and  1.229 
for  y  -  1.1058.  1 

This  solution,  or  any  other  such  simple  solution,  is  not  adequate  here 
because 

(1)  The  gap  between  the  pipes  is  relatively  narrow  in  comparison  with  the 
pipe  dimensions,  and  backpressure  on  the  outside  of  the  process  pipe  must  be 
accounted  for  in  a  force  balance; 

(2)  The  non-uniform  filling  of  the  gap  between  the  process  and  guard 
pipes  with  steam  and  compressed  air  makes  ill-defined; 

(3)  The  nozzle  discharge  coefficient  is  not  known  a  priori  to  obtain  the 
effective  value  of  Ag  for  the  break  in  the  process  pipe; 

(4)  Time-dependent  pressures  must  be  considered  here  initially  to  obtain 
the  early-on  transient  airshock/steam  shock  contributions  to  the  guard  pipe 
loading; 

*Moody,  F,  J.,  "Prediction  of  Blowdown  Thrust  and  Jet  Forces,"  ASME  Paper  69-HT-31, 
ASME-AIChE  Heat  Transfer  Conference,  Minneapolis,  Minn.,  Aug  1969. 
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(5)  The  distribution  of  pressure  along  the  guard  pipe  inside  surface, 
not  just  the  total  force,  is  needed  for  future  structural  response  calculations. 

The  total  force  calculations  are  useful,  however,  for  comparing  with 
integrated  results  obtained  from  the  hydrocode  computations.  These  comparisons 
are  discussed  later. 
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CHAPTER  3 

HYDROCODE  SOLUTIONS 


This  was  a  problem  well  suited  for  an  Eulerian  (rather  than  Lagrangean) 
hydrocode  since  there  is  much  distortion  in  the  flow.  In  the  initial  phase 
of  the  analysis  two  hydrocodes  were  used,  TUULI  and  CSQ. 

TUULI*  is  a  two-dimensional  Eulerian  hydrocode,  written  in  1975  at  NSWC 
by  D.  Lehto.  Full  documentation  is  not  yet  available.  The  code  is  based  on 
the  fluid-in-cell  (FLIC)  method  (REFERENCE  2)  where  the  calculation  is  done 
in  four  steps  per  time  cycle: 

1)  The  accelerations  are  calculated  from  the  pressure  gradients  and  new 
provisional  velocities  are  calculated  (without  convection). 

2)  The  provisional  velocities  are  used  to  calculate  the  pdV  work  done  on 
each  zone;  this  gives  provisional  internal  energies  (still  without  convection). 

3)  The  material  transport  (convection)  is  done  with  the  provisional 
velocities  and  energies. 

4)  Any  adjustments  ..jeded  to  conserve  both  energy  and  momentum  are  made. 
Any  kinetic  energy  correction  needed  for  momentum  conservation  is  taken  from 
(or  added  to)  the  internal  energy.  These  adjustments  are  necessary  because 
the  flow  mixes  dissimilar  flows  from  adjacent  zones,  and  both  kinetic  energy 
and  momentum  cannot  be  conserved  in  a  mixing  process.  This  fourth  step  is 
explicitly  done  in  TUULI  because  internal  energy  is  transported;  it  is 
implicitly  done  in  the  original  FLIC  code  because  total  energy  is  transported. 
This  is  a  trivial  arbitrary  choice. 

TUULI  handles  two  materials  (here,  steam  and  air)  and  has  an  option  for 
assigned  inflow  in  any  chosen  zones  (used  extensively  in  these  calculations). 

The  grid  is  composed  of  fixed  rectangular  zones.  Shocks  are  handled  by  the 
quadratic  artificial  viscosity  method. 

- - 

Gentry,  R.  A,,  Martin,  R.  E.,  and  Daly,  E.  J.,  "An  Eulerian  Differencing 
Method  for  Unsteady  Compressible  Flow  Problems,"  J.  Comput.  Phys.  1_,  pp  87-118, 
1966. 

♦Until  now,  this  code  was  called  TUTTI,  However,  this  is  also  the  name  of  a  Los 
Alamos  Scientific  Laboratory  equation-of-state  program  (J.  Appl.  Phys.  5^  (10) 
5360  (1980)).  Thus,  the  name  change:  TUTTI  to  TUULI. 
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CSQ  is  a  two-dimensional  Eulerian  hydrocode  written  at  the  Sandia  Laboratories. 
Documentation  for  this  code  is  provided  in  Reference  3.  In  this  method,  the 
calculation  is  performed  in  two  main  phases: 

1)  The  finite  difference  analogs  of  the  complete  two-dimensional  Lagrangean 
equations  are  solved  during  each  time  cycle. 

2)  At  the  end  of  the  time  cycle,  the  code  rezones  the  mesh  back  to  the 
original  configuration. 

The  net  result  of  the  rezoning  is  an  Eulerian  calculation.  The  version  of 
CSQ  described  in  Reference  3  handles  only  two  materials  in  any  chosen  zone. 

The  grid  is  composed  of  fixed  rectangular  zones.  Shocks  are  handled  by  the 
quadratic  artificial  viscosity  method. 

The  computations  were  performed  on  a  CYBER  176  at  the  Air  Force  Weapons 
Laboratory  in  Albuquerque,  New  Mexico  over  telephone  lines  from  NSWC. 

3.1  ZONE  SIZE  SELECTION 

Both  TUULI  and  CSQ  face  the  same  limitations  on  available  computer  size 

and  cost  of  computer  time.  Care  was  taken  to  choose  a  calculation  mesh  just 

fine  enough  to  give  the  desired  accuracy;  this  was  done  by  simply  trying 

progressively  finer  meshes  until  the  overpressure  loading  on  the  guard 

pipe  inside  wall  no  longer  changed  significantly.  The  final  mesh  selected 

was  0.315  zones/mm  (8  zones/in)  within  the  channel  formed  by  the  break  after 

it  was  determined  that  reduction  of  the  mesh  size  down  to  0.630  zones/mm 

(16  zones/in)  produced  no  significant  effect  on  the  calculated  flow. 

The  mesh  was  nonuniform  and  varied  for  each  problem  geometry,  being  fine  in 

the  channel  and  near  the  impact  area  on  the  guard  pipe  and  progressively 

coarser  with  increasing  distance  from  the  region  of  interest. 

» 

3.2  ZONE  SHAPE  LIMITATION 

Both  hydrocodes  are  limited  to  rectangular  zones.  This  makes  the 
longitudinal-break  problem  awkward,  because  the  concentric  circles  needed 
to  represent  the  problem  cross  section  would  have  to  be  made  up  of  rectangular 
steps  (See  Figure  2).  For  the  longitudinal -break  calculations,  the  pipes  were 
straightened  out  into  parallel  planes  as  shown  in  the  figure;  this  is  expected 
to  be  a  good  approximation  since  the  separation  distances  between  pipes  are 
relatively  small  compared  with  the  pipe  radii  for  the  geometries  investigated 
here.  No  problem  arose  in  setting  up  the  mesh  with  rectangular  zones  for  the 
circumferential -break  calculations  because  all  bounding  surfaces  could  be 
represented  by  straight  lines.  The  computing  mesh  for  the  flat-plate 
approximation  for  the  longitudinal  break  shown  in  Figure  2  looks  similar  to 
the  mesh  for  the  circumferential  break;  however,  the  latter  has  axial  symmetry. 


3 

Thompson,  S.  L.,  "CSQ  --  A  Two  Dimensional  Hydrodynamic  Program  with  Energy 
Flow  and  Material  Strength,"  Sandia  Labs.  SAND  74-0122,  Aug  1975. 
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3.3  EQUATIONS  OF  STATE 

When  the  steam  pressure  drops  as  it  leaves  the  ruptured  process  pipe, 
condensation  takes  place  and  the  resulting  two-phase  flow  expands  isentropically. 
For  example,  the  isentropic  exponent  for  equilibrium  flow  (i.e.,  equilibrium 
between  the  vapor  and  droplet  phases)  is  1.1058  (Table  2)  for  saturated  steam  at 

8.3  MPa  (1200  psia).*  However,  it  takes  a  finite  time  for  the  vapor  to  condense 
into  droplets.  If  this  relaxation  were  a  relatively  slow  process  with  respect 
to  the  expansion,  the  expanding  steam  would  continue  to  behave  as  a  pure  vapor. 

On  the  vapor  side  of  the  saturation  line,  the  adiabatic  exponent  is  1,2592; 

a  frozen  flow  would  be  expected  tn  maintain  this  exponent  beyond  the  saturation 
line.  The  actual  flow  would  lie  between  these  two  states  during  the  expansion 
process.  Referring  to  Figure  3  (taken  from  Reference  5),  the  Tl-curve  is 
the  equilibrium  flow  isentrope  and  the  I '-curve  segment  below  the  saturated 
vapor  curve  (SV)  is  the  frozen  flow  isentrope.  The  actual  curve,  nonequilibrium 
flow  designated  ACT,  for  the  two-phase  steam  mixture  which  takes  into  account 
the  kinetics  of  the  relaxation/condensatiori  process  falls  somewhere  in  between 
the  TI  and  the  I1  curves  as  shown  in  the  figure. 

The  effect  of  nonequilibrium  flow  (delay  in  condensation)  was  bounded 
by  performing  hydrocode  calculations  with  adiabatic  exponents  {y)  for  both 
the  frozen  and  the  equilibrium  flow  conditions.  An  ideal-gas  equation  of 
state. 


P  •  (y  -  1)pE  , 

was  used  to  give  the  isentrope, 

Pp"Y  =  constant, 

where  P  is  pressure,  p  is  density,  and  E  is  internal  energy.  The  results, 
given  in  Figure  4  for  a  typical  circumferential-break  geometry,  indicate 
that  nonequilibrium  flow  corrections  to  the  pressure  loading  on  the  guard 
pipe  opposite  the  pipe  break  are  negligible.  Equilibrium  flow  was  assumed 
for  the  rest  of  the  hydrocode  calculations. 


♦Steam  table  (Reference  4)  pressure-volume  data  were  fitted  (on  a  log-log 
scale)  with  straight  lines  to  get  the  adiabatic  exponents  given  in  Table  2. 

^Keenan,  0.  H.  and  Keyes,  F.  G.,  Thermodynamic  Properties  of  Steam.  Including 
Data  for  the  Liquid  and  Solid  Phases,  1st' ed. ,  Wiley,' New  York,  1936. 

5 

Zel'dovich,  Ya.  B.,  and  Raizer,  Yu,  P.,  Physics  of  Shock  Waves  and  High 
Temperature  Hydrodynamic  Phenomena,  Volume  II,  ed.  by  U.  b.  Hayes  and 
R.  F.  Probstein,  Academic  Press,  New  York,  1967. 
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TABLE  2 

saturated  steam  uata 


METRIC  UNITS- 


PMESSURE 

(MPA) 

b.b*S 

b  •  2  74 

b.blb 

SPECIFIC  VOLUME 

(M3/MG) 

27. H2 

22.6V 

21.54 

i)EnS I  fY 

(MG/M  3) 

O.OJb*b 

0.04426 

0.04644 

TEMPEkATURE 

(K) 

5S/.9B 

S70.S1 

573.39 

tNTROPY 

( J/KG.K) 

3blb. 

S722. 

5699. 

ADIABATIC  CONSTANT 

ON  VAPOR  SIDE 

1.2623 

1.25*2 

1.2777 

ADIABATIC  CONSTANT 
ON  2-PHASE  SIDE 

1  •  0866 

1.1056 

i .104b 

t  NGL  I  SH  UNITS  — 
KMESSURE 

(PS  I  A ) 

1000 

1200 

1230 

SPEC  IK IC  VOLUMt 

(FT  3/LH ) 

0.44bb 

0.3b  IV 

0. 3460 

Tt^PtHATUPE 

( F  > 

344.67 

S67.22 

5  7  2  •  4  2 

ENTROPY 

(BTU/F .LB) 

1 .3H*7 

1.3667 

1 .3bl2 

T ABLE  3 

tAlT  CONDITIONS  AT  SLIT 


b.274  MPA  <1200  PblA)  SATURATED  5  T  t AM ,  G AMM A  =  1.1054 


bkeak  wioth« 

b0.«  MM  (? 

IN) .  PIPE 

thick n tbs 

=50. b  MM 

(2  IN)  . 

distance  FROM 

pakt icle 

5UOND 

MASS 

MluPLANE  OF 

PKtSSURE 

UtNsi TY 

VELOCITY 

sPtED 

FLUX 

BhEAK  (MM) 

(MPA) 

( G/MM3 ) 

(M/S) 

(M/S) 

(G/M2/US) 

1.6 

2.613 

lb,  52 

6UV.H 

431.5 

*.464 

4.8 

2.616 

15.32 

hOb.5 

431 . 7 

*  .444 

/.9 

2.61b 

15.32 

605.  V 

431  .  V 

*.404 

11.1 

2.617 

15.52 

601.1 

431  .B 

*.32* 

14.3 

2.606 

lb. 52 

3*1.7 

430.  V 

*.  lbj 

17.  b 

2.S73 

15.52 

573.2 

426.2 

b.b*  6 

20.6 

2.506 

15.52 

335, H 

422.6 

a.  316 

23. B 

2.403 

15.52 

463.7 

4l3.b 

7.19  7 

SONIC 

27.6b 

472. B 

*♦72.6 

12.05 

16 


FIGURE  4.  EFFECT  OF  THERMODYNAMIC  NONEQUILIBRIUM  ON  PRESSURE  ON  GUARD  PIPE 
OPPOSITE  THE  BREAK 
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The  air  in  the  gap  between  the  process  and  guard  pipes  was  treated  as 
an  ideal  gas  with  y  =  1.4.  The  air  in  the  gap  was  initially  set  to  an  average 
temperature  of  422  K  (300°  F)  for  all  the  hydrocode  calculations.  This  value 
was  determined  by  performing  a  steady  one-dimensional  heat  conduction  analysis 
(see  Reference  6,  p.  37)  for  a  steam  pipe  surrounded  by  an  air  gap,  a  guard 
pipe,  and  an  ambient  atmosphere  with  free  convection. 

3.4  HYDROCODE  SELECTION 


TUULI  and  CSQ  give  essentially  the  same  results,  as  expected.  This  is 
demonstrated  in  Figure  5  which  presents  the  pressure  versus  time  history  from 
each  code  at  the  guard  pipe  surface  directly  opposite  the  circumferential  break 
for  Run  No.  4  listed  in  Table  1.  The  small  oscillations  are  unresolved  airshock 
reflections  between  the  nose  of  the  steam  jet  and  the  guard  pipe  wall.  They 
are  considered  insignificant  in  terms  of  pipe  response  to  the  pressure  loads. 


TUULI  was  the  hydrocode 
because  the  authors  are  more 


chosen  for  the  remainder  of  the  calculations 
familiar  with  it  than  with  CSQ. 


^Kreith,  F. ,  Principles  of  Heat.  Transfer,  2nd  ed.,  International  Textbook  Co., 
Scranton,  Pa.,  1965. 
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FIGURE  5.  TUULI  AND  CSQ  CODE  PRESSURE  VS  TIME  ON  GUARD  PIPE  OPPOSITE  THE  BREAK 


.1 


NSWC  TP  80-229 


CHAPTER  4 
HYDROCODE  RESULTS 


The  expected  flow  phenomena  following  an  instantaneous  break  in  the  process 
pipe  (for  the  ci rcumferential -break  geometry)  are  shown  in  Figure  6.  A 
rarefaction  wave  recedes  into  the  process  pipe  steam  reservoir.  A  steam  plume 
jets  from  the  break  in  the  process  pipe  into  the  air  gap  between  the  process 
and  guard  pipes.  Sonic  conditions  are  attained  within  the  break  and  the  flow 
becomes  choked.  An  airshock  precedes  the  steam  jet  and  a  stationary  shock 
front  (in  the  steam)  is  formed  which  stands  off  at  some  distance  from  the 
guard  pipe  inner  wal 1 . 

In  section  4.1  several  elements  of  the  flow  phenomena  --  rarefaction  wave, 
sonic  flow  in  the  process  pipe  break,  and  discharge  coefficient  for  the  break  -- 
are  discussed.  These  elements  allow  a  substantial  simplification  of  the  flow 
field  calculations.  The  flow  internal  to  the  process  pipe  need  only  be  computed 
once  (with  a  hydrocode)  for  each  initial  steam  flow  condition.  Following  this, 
truncated  flow  field  computations  are  then  performed  for  the  different  process/ 
guard  pipe  geometries.  The  truncated  flow  field  calculation  involves  specifying 
the  steady  outflow  from  the  process  pipe  exit  and  then  computing  the  flow  field 
for  the  intervening  space  between  the  process  pipe  and  the  guard  pipe. 

Section  4.2  discusses  the  following  aspects  of  the  flow  field  for  the 
region  between  the  process  pipe  and  the  guard  pipe: 

1)  The  pressure  loads  on  the  guard  pipe  wall  produced  by  the  airshock 
preceding  the  steam  plume; 

2)  The  transient  initial  pressure  loads  produced  by  the  steam  plume  on 
the  guard  pipe  wal 1 ; 

3)  The  stationary  shock  in  the  vicinity  of  the  guard  pipe  wall;  and 

4)  The  maximum  pressure  loads  on  the  guard  pipe  wall  that  are  established 
by  the  steady  flow  of  the  steam  exiting  from  the  process  pipe  break. 

4.1  OUTFLOW  FROM  PROCESS  PIPE 

RAREFACTION  WAVE.  A  rarefaction  wave  runs  from  the  break  into  the  steam  in 
the  process  pipe.  For  a  circumferential  break,  this  wave  reflects  from  the  pipe 
axis  and  part  of  it  comes  back  into  the  break.  Figure  7  shows  the  guard  pipe 
wall  reflected  pressure  versus  time  for  such  a  calculation  carried  out  far 
enough  in  time  to  follow  the  propagation  of  the  reflected  rarefaction  wave  to 
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FIGURE  7.  EFFECT  OF  AXIAL  REFLECTION  OF  RAREFACTION  WAVE  IN  PROCESS  PIPE 
ON  GUARD  PIPE  PRESSURE 
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and  into  the  break.  The  net  effect  on  the  guard  pipe  reflected  pressure  loading 
is  a  negligible  pressure  oscillation  that  arrives  after  an  essentially  steady 
flow  pattern  at  the  guard  pipe  wall  has  been  attained.  The  shape  of  the 
pressure  oscillation  indicated  in  Figure  7  represents  the  interaction  of  the 
rarefaction  wave  with  the  stationary  shock  and  the  guard  pipe  wall  --  the  lack 
of  definition  of  this  pressure  oscillation  is  caused  by  the  coarse  zone  size 
used  in  this  computation. 

SONIC  FLOW  IN  PROCESS  PIPE  BREAK.  Within  the  break,  sonic  flow  conditions 
occur  at  the  minimum  flow  cross  section  if  the  downstream  pressure,  p^,  is 


where  pu  is  the  upstream  (reservoir)  pressure  and  y  is  the  specific  heat  ratio 
of  the  gas. 

For  a  free  jet,  the  downstream  pressure  is  equal  to  the  ambient  pressure. 

For  a  confined  jet  such  as  discussed  here,  the  downstream  pressure  is  not 
defined  a  priori  because  of  the  pressure  buildup  in  the  confined  space.  This 
does  not  present  a  problem  for  the  hydrocode  calculation  because  it  automatically 
simulates  the  correct  flow,  whether  it  is  sonic  or  subsonic.  However,  if 
preliminary  hydrocode  calculations  indicate  that  the  flow  in  the  break  remains 
sonic  for  the  duration  under  investigation  (which  turns  out  to  be  the  situation 
here),  then  considerable  savings  in  computer  costs  can  be  made  by  assigning 
the  proper  steady  flow  conditions  at  the  break  exit  and  calculating  only  the 
flow  downstream  of  the  break  for  the  various  pipe  spacings.  The  downstream 
flow  cannot  send  any  signals  through  the  sonic  f 1 ov/  into  the  process  pipe. 

DISCHARGE  COEFFICIENT  FOR  PROCESS  PIPE  BREAK.  The  process  pipe  break  is 
expected  to  have  sharp  edges.  For "such"  a  nozzle,"  the  flow  contracts  and 
does  not  fill  the  "nozzle"  cross  section  (Figure  8).  This  effect  is  the 
familiar  vena  contracta  of  incompressible  flow  and  is  known  tc  occur  for 
compressible  flow  as  well  (References  7  -  11).  The  area  contraction 

^Bean,  H.  S.,  Buckingham,  E.  and  Murphy,  P.  S.,  "Discharge  Coefficients  of 
Square-Edged  Orifices  for  Measuring  the  Flow  of  Air,"  N.B.S.J.  Res.  2, 
pp  561-568,  1969. 

^Stanton,  T.  E.,  "On  the  Flow  of  Gases  at  High  Speeds,"  Proc.  Roy.  Soc.  Ill , 
pp  306-339,  1926. 

^Perry,  J.  A.,  Or.,  "Critical  Flow  Through  Sharp-Edged  Orifices,"  Trans.  Am. 

Soc.  Hech.  Engrs.  71 ,  pp  757-764,  1949. 

^Grace,  H.  P.  and  Lapple,  C.  E. ,  "Discharge  Coefficients  of  Small -Pi ameter 
Orifices  and  Flow  Nozzles,"  Trans.  ASME,  pp  639-647,  Jul  1951. 

^Arnberg,  B.  T. ,  "Review  of  Critical  Flowmeters  for  Gas  Flow  Measurements," 

J.  Basic  Engq. ,  pp  447-46'',  .tec  1362. 
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coefficient  is  typically  between  0.6  and  1.0.  This  contraction  allows  the 
downstream  pressure  to  have  some  influence  on  the  outflow  even  though  the  flow 
is  sonic,  because  signals  can  reach  the  sonic  region  by  bypassing  the  sonic 
core  of  the  flow.  Because  no  data  were  available  on  the  flow  of  steam  through 
sharp-edged  channels,  the  "nozzle"  or  process  pipe  break  flow  was  determined 
as  part  of  the  problem.  The  hydrocode  calculations  gave  a  discharge 
coefficient  of  0.77,  compared  with  1.00  expected  for  a  rounded-inlet  orifice. 
This  discharge  coefficient  is  used  in  calculating  the  outflow  from  the  process 
pipe  break  exit  for  some  of  the  truncated  flow  field  calculations. 

The  steam  jet  reaches  the  end  of  the  50.8  mm  (2.00  in)  long  break  (process 
pipe  wall  thickness)  before  the  flow  has  expanded  to  fill  the  entire  cross- 
sectional  area  for  the  50.8  mm  break  width.  This  result  is  evident  in  Table  3 
which  gives  the  steady-state  outflow  parameters,  computed  using  CSQ,  for  the 
8.3  MPa  (1200  psia)  process  pipe  reservoir  conditions.  The  mass  flux  near  the 
centerline  is  30%  greater  than  the  mass  flux  near  the  "nozzle"  wall.  The 
effect  is  also  indicated  schematically  in  Figure  8  for  the  "nozzle"  with  the 
sharp-edged  inlet  --  the  figure  indicates  that  the  supersonic  flow  has  not 
expanded  to  encompass  the  entire  area  of  the  "nozzle,"  Results  such  as  given 
in  Table  3  will  be  used  to  avoid  calculation  in  every  computer  run  of  the 
interior  of  the  process  pipe  and  the  interior  of  the  break. 

TRUNCATED  FLOW  FIELD  CALCULATIONS.  To  reduce  computer  time  costs,  the 
runs  with  8.3  MPa  (1200  psia)  process  pipe  reservoir  conditions  (Table  1, 

Runs  No.  1  -  12,  16,  17),  were  computed  using  a  truncated  flow  field.  The 
calculation  of  the  flow  in  the  interior  cf  the  process  pipe  (the  reservoir) 
and  the  interior  of  the  break  (the  "nozzle")  was  replaced  by  the  steady-state 
flow  conditions  of  Table  3  at  the  break  outlet  since  the  effect  of  the  rare¬ 
faction  wave  moving  into  the  reservoir  was  shown  to  be  negligible  in  Figure  7. 

For  the  12.7  mm  (0.5  in)  wide  breaks  in  these  runs,  the  flow  filled  the 
break  exit,  so  the  exit  conditions  were  obtained  from  analytical  sonic  flow 
calculations  over  0.77  (the  discharge  coefficient  determined  for  these  sharp- 
edged  breaks  or  "nozzles")  of  the  break  area  expanded  to  the  full  exit  area 
with  the  usual  shock-tube  area  relations.  The  sonic  conditions  given  in  the 
last  line  of  Table  3  apply  across  the  entire  exit  plane  of  the  12.7  mm  (0.5  in) 
wide  breaks. 

4.2  PRESSURES  ON  GUARD  PIPE  WALL 

AIRSHOCK  FORMATION  BETWEEN  PIPES.  The  exiting  steam  jet  drives  the  air 
between  the  pipes  ahead  of  it,  and  forms  at  first  a  bow  airshock  at  the  nose 
of  the  steam  plume,  and  later  an  airshock  that  propagates  ahead  of  the  steam 
jet  flow  in  the  space  between  the  pipes. 

Figure  9  presents  a  typical  pressure  distribution  along  the  guard  pipe 
inside  wall  showing  the  presence  of  the  airshock.  The  airshock  shown  in 
the  figure  has  an  overpressure  of  0.5  to  1.0  MPa  (5  to  10  bars)  and  is  not 
regarded  as  structurally  significant;  it  is  therefore  handled  with  coarse 
zoning. 
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FIGURE  9.  PRESSURE  DISTRIBUTION  ALONG  GUARD  PIPE,  INCLUDING  AIR  SHOCK 
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In  the  remainder  of  this  report,  the  discussion  concentrates  on  the  sharp 
pressure  maximum  on  the  guard  pipe  inside  surface  facing  the  site  of  the  process 
pipe  break. 

PRESSURE  TRANSIENTS  ON  THE  GUARD  PIPE  HALL.  Before  discussing  the  main  results, 
the  pressure  transients  due  to  shock  reflections  between  the  oncoming  steam  jet 
and  the  guard  pipe  are  examined  briefly.  These  are  poorly  resolved  in  the 
two-dimensional  calculations  (e.g.,  see  the  oscillations  in  Figure  5).  To 
resolve  these  shocks,  the  zoning  scale  would  have  to  be  about  an  order  of 
magnitude  finer  and  the  computer  cost  for  a  two-dimensional  calculation  would 
be  quite  high.  However,  an  upper  bound  can  be  obtained  with  a  finely-zoned 
one-dimensional  (WUNDY  code;  Reference  12)  Lagrangean  calculation.  The 
simple  linear  shock  tube  geometry  and  the  relationship  with  the  geometry  of 
the  real  problem  are  shown  in  Figure  10.  The  resulting  pressure-time  record 
at  the  guard  pipe  is  shown  in  Figure  11.  The  reflected  shocks  are  not 
significant  because: 

1)  They  are  very  short  (about  10  ps)  and  carry  litt'e  impulse; 

2)  They  would  be  reduced  by  the  two-dimensional  expansion  for  the  real 
problem  geometry; 

3)  They  would  be  very  much  reduced  in  amplitude  by  the  finite  opening 
time  of  the  process  pipe  break  which  is  assumed  to  be  instantaneous  for  these 
calculations. 

The  opening  time  of  a  circumferential  50.8  nm  (2.00  in)  break  in  an 
infinite  pipe  with  the  parameters  of  Runs  1  through  12  is  about  27  ms  (Appendix 
A);  the  opening  time  of  a  longitudinal  break  is  about  8  ms,  calculated  using 
a  clamshell -type  opening  model  with  hinging  at  the  side  opposite  the  break. 

These  long  opening  times  mean  that  the  steam  jet  that  first  reaches  the 
guard  pipe  originated  from  a  small  slit  and  would  be  greatly  attenuated  in 
pressure  because  the  distance  to  the  guard  pipe  would  represent  many  slit 
widths. 

These  insignificant  transient  effects  were  ignored.  The  computational 
effort  was  directed  towards  determining  the  steady-state  pressure  distribution 
on  tne  guard  pipe.  The  slit  (process  pipe  break  exit)  was  assumed  to  open 
instantaneously,  which  allowed  rapid  attainment  of  the  same  steady  flow  that 
would  have  been  eventually  reached  with  a  slowly-opening  slit. 

The  steady-flow  pressures  reported  in  the  next  section  represent  over¬ 
estimates  if  the  pipes  are  so  short  that  the  supply  pressure  drops  before 
steady  flow  is  attained. 
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FIGURE  10».  GEOMETRY  OF  REAL  PROBLEM 
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FIGURE  10b.  GEOMETRY  OF  LINEAR  SHOCK  TU3E  PROBLEM 
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STEADY-FLOW  PRESSURES  ON  GUARD  PIPE  WALL.  The  main  results  of  this  study 
are  presented  In  Figures  12  -  lf>.  These  figures  present  the  steady-flow 
pressure  distribution  on  the  guard  pipe  inside  surface  in  the  vicinity  of  the 
process  pipe  break.  The  full  distribution  of  computed  pressure  data  including 
the  airshock  effects  are  listed  in  corresponding  Tables  4-6.  The  longitudinal 
break  maximum  pressures  are  greater  than  the  corresponding  circumferential -break 
maximum  pressures,  because  the  latter  flow  has  more  geometric  spreading  (compare 
Figures  12  and  13  and  see  Figure  14). 

In  both  process  pipe  break  configurations,  the  maximum  pressure  drops  off 
rapidly  with  increasing  process/guard  pipe  separation  distances. 

A  check  was  made  on  how  the  total  force  on  the  guard  pipe  plus  the  process 
pipe  agrees  with  the  thrust  equation  for  a  jet  of  outlet  area  A  hitting  an 
unconfined  plate  (Reference  1).  Repeating  the  expression  given  in  Chapter  2, 

T/Ae  -  ClP0  -  P» 

with  «  1.229  (for  y  »  1.1056),  P0  =  8.274  MPa  and  P*  =  0.101  MPa 
gives  T/Afi  -  1.017  x  108  dyne/cm2  *  10.17  M/m2. 

This  can  be  compared  with  the  net  outward  force  on  both  pipes  for  the 
parallel-plate  model  of  the  longitudinal  break: 


Fout‘/Vfl  -  /Vdfl 

gp  Pp 

For  small  plate  separations,  this  integration  removes  most  of  the  effect  of 
the  shock  running  between  the  plates,  because  the  shock  pressure  distributions 
are  nearly  identical  on  the  two  surfaces.  For  large  plate  separations,  the 
shock  contribution  does  not  cancel. 

Figure  16  shows  that  the  calculated  net  forces  are  close  to  those  of  the 
thrust  equation.  These  integrals  should  not  be  regarded  as  precise,  because 
they  depend  strongly  on  results  in  the  rather  coarsely  zoned  airshock  region. 
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FIGURE  12 .  STEADY-FLOW  PRESSURE  ALONG  GUARD  PIPE  FOR  50.8  mm  (2  in) 
CIRCUMFERENTIAL  BREAK 


NSWC  TR  80  229 


i 


i 


FIGURE  13.  STEADY-FLOW  PRESSURE  ALONG  GUARD  PIPE  FOR  60.8  mm  (2  in) 
LONGITUDINAL  BREAK 
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TABLE  4 

STEADY  FLOW  PRESSURES  ALONG  INSIDE  Of  GUARD  PIPE  DUE  TO  STEAM  JET 
FROM  50.8  MM  (2  IN)  RREAK  IN  PROCESS  PIPE 
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Z  =  UISTmNCE  ALUNO  GUAM!)  PIPE  f  HUH  CENlTfcR  OF  bRtAK  (MM). 
F=AbSOuUTE  PRESSURE  On  bUAPU  HlPt  (MPA) . 
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TAHLt  b 

STEADY  FLOW  PRESSURES  ALONG  INSIDE  OF  GUARD  PIPE  DOE  TO  STEAM  JET 
FROM  12.7  MM  (0.5  IN)  BREAK  11!  PROCESS’ PIPE 
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TAHLfc  6 

STEADY  FLOW  PRESSURES  ALONG  INSIDE  OF  GUARD  PIPE  DUE  TO  STEAM  JET 

FROM  BREAK  IN  PROCESS  PIPE 
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Figure  17  shows  the  pressure  distributions  (indicating  the  presence  of 
the  stationary  shock)  along  the  plane  of  symmetry  passing  through  the  center 
of  the  longitudinal  break  for  selected  process/guard  pipe  separation  distances. 
The  stationary  shock  (in  the  steam)  shows  the  characteristic  smoothing-out  of 
the  discontinuity  that  is  inherent  in  the  artificial-viscosity  method  for 
handling  shock  discontinuities  numerically. 

Some  details  cf  the  steady  flow  (for  longitudinal-break  configuration) 

In  Runs  7  and  11  are  presented  in  Figures  18  and  19,  respectively.  These 
runs  differ  only  in  the  spacing  between  the  pipes  (4  and  6  inches,  respectively) 
The  solid  curves  are  isobars.  There  is  a  region  of  backflow  towards  the  break 
in  both  figures.  The  entire  flow  field  in  Figure  18  is  composed  of  steam. 

Figure  19  shows  a  trapped  streamer  of  air,  which  conveniently  delineates  the 
boundary  of  the  steam  jet. 
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FIGURE  18.  FLOW  PATTERN  FOR  PROBLEM  NO.  7874  (RUN  7) 
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FIGURE  19 .  FLOW  PATTERN  FOR  PROBLEM  NO.  7876  (RUN  11) 
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CHAPTER  5 
CONCLUSIONS 


The  results  in  Figures  12  -  15  indicate  that  the  steady-flow  maximum 
pressure  loads  on  the  guard  pipe  inside  surface  produced  by  ruptures 
(longitudinal  and  circumferential)  of  the  steam  process  pipe  are  less  than 
the  Initial  process  pipe  steam  pressure  (as  Bernollis1  law  would  indicate). 

As  the  process/guard  pipe  separation  distance  approaches  zero  the  pressure 
loading  approaches  the  initial  process  pipe  steam  pressure.  The  guard  pipe 
pressure  loads  fall  off  rapidly  with  distance  from  the  process  pipe  break 
site  along  the  guard  pipe  surface.  Transient  airshocks  with  high  peak 
pressures  (and  very  short  durations)  occur  before  the  flow  becomes  steady; 
however,  these  airshocks  do  not  transmit  sufficient  momentum  to  the  guard 
pipe  wall  to  affect  the  guard  pipe  structural  response  appreciably.  A  guard 
pipe  designed  to  withstand  the  process  pipe  pressure  can  withstand  a  steam 
jet  from  a  break  in  the  process  pipe.  The  effect  of  whipping  of  the  process 
pipe  is  not  considered  in  this  report. 

Only  simple  size  scaling  can  be  accurately  applied  to  these  results 
(i.e.,  multiplying  all  dimensions  and  time  by  the  same  factor).  The  maximum 
pressure  on  the  guard  pipe  depends  strongly  on  the  ratio  of  distance  to  guard 
pipe  divided  by  break  width,  and  weakly  on  the  ratio  cf  break  length  (i.e., 
pipe  thickness)  to  break  width  and  on  the  equation  of  state  for  steam. 
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APPENDIX  A 

OPENING  TIME  OF  CIRCUMFERENTIAL  BREAK  IN  PROCESS  PIPE 


Assume  that  the  break  occurs  instantaneously  and  calculate  how  long  it 
takes  for  the  ends  to  pull  apart  to  the  desired  break  size.  The  mechanism 
for  pulling  the  ends  apart  is  a  relief  wave  that  propagates  along  the  pipe 
and  drops  the  meridional  stress  from  its  loaded  condition  to  zero  stress. 
The  meridional  stress  is 


0  '  =  36-20  MPa  s  5250  psi 

where  p  =  pressure  in  pipe  (=8.274  MPa  =  1200  psia) 
R  =  radius  of  pipe  (=0.445m  =  17.5  in) 
t  3  thickness  of  pipe  (=50.0  mm  =  2  in) 


The  speed  of  a  longitudinal  wave  in  an  infinite  plate  is  (Reference  A1  -  p.  81) 


L  fp (1- 


T  =  5.24  x  103  m/s 
vZ)  - 


where  C  =  Young's  modulus  (=2.0  x  105  MPa  for  steel  (Reference  A2  -  p.  2-68)) 

v  =  Poisson's  ratio  (=0.3) 

p  =  density  (=8  x  103  kg/m3) 


The  speed  of  the  pipe  (particle  velocity)  in  the  stress-relieved  region  is 


^Ko'lsky ,  H.  > 

A2Gray,  D.  E. 
1972. 


v 


0.948  m/s 


Stress  Waves  in  Solids,  Dover,  New  York,  1963. 

(ed.),  American  Institute  of  Physics  Handbook,  3rd  ed.,  McGraw-Hill, 
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To  open  a  50.8  min  (2  in)  break,  each  end  has  to  move  25.4  mm  (1  in).  The 
time  this  takes  is 

t  (25.4  mm  s  1  in)  •  (25.4  mm)  /  (0.948  x  103  mm/s)  ■  26.8  ms 


In  this  time,  the  elastic  relief  wave  has  moved  down  the  pipe  a  distance  of 
(0.0268  s)  x  (5.24  x  105  m/s)  ■  140  m.  Anything  within  70  m  that  inhibits 
longitudinal  motion  of  the  pipe  would  send  back  a  wave  that  would  slow  the 
opening  of  the  gap  before  it  reached  a  total  size  of  50.8  mm  (2  in). 
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